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The unresolved hard X-ray background: the missing source 
population implied by the Chandra and XMM-Newton deep fields 
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ABSTRACT 

We extend our earlier work on X-ray source stacking in the deep XMM-Newton observation 
of the Lockman Hole, to the 2 Ms Chandra Deep Field North and the 1 Ms Chandra Deep 
Field South. The XMM-Newton work showed the resolved fraction of the X-ray background 
to be ~ 80-100 per cent at < 2 kcV but this decreased to only ^ 50 per cent above ~ 8 kcV. 
The CDF-N and CDF-S probe deeper, and are able to fill-in some of the missing fraction in 
the 4-6 keV range, but the resolved fraction in the 6-8 keV band remains only ^ 60 per cent, 
confirming the trend seen with XMM-Newton. The missing X-ray background component has 
a spectral shape that is consistent with a population of highly obscured AGN at redshifts 
^ 0.5-1.5 and with absorption column densities of ^ lO^'^-lO^'' cm^^. 

Key words: surveys - galaxies: active - cosmology: diffuse radiation - X-rays: galaxies - 
X-rays: diffuse background 



1 INTRODUCTION 

The nature of the X-ray background (XRB) has been a subject of 
contention since its discovery 40 years ago (Giacconi et al. 1963). 
It is now clear that that at energies above ^ 1 keV the background 
is made up of the summed emission from point sources, princi- 
pally Active Gal actic Nuclei (AGN). T he spectrum measured by the 
HEAO satellite jMarshall et alll98(t) is similar in shape to 40 keV 
bremsstrahlung and is well-described by a F — 1.4 power-law from 
~ 1 up to ~ 15 keV where it starts to flatten off before peaking in 
vli, at ~ 30 keV. Below ~ 1 keV the background is dominated by 
diffuse Galactic and local bubble emission; the purely extragalactic 
component has not been well determined. 

Resolving the XRB into discrete sources is crucial to under- 
standing the nature of the AGN and galaxy populations in the 
Universe. It also allows us to place limits on the residual XRB 
emission which can be attributed to other processes, including that 
from truly diffuse extragalactic emission. A significant amount 
(~ 70-80 per cent) of the 0.5-2 keV XRB was resolved into point 
sources by ROSAT ( Hasinaer et al. 1998) but progress at harder en- 
ergies has required the improved sensitivity and spatial resolution 
of the Chandra and XMM-Newton observatories. XMM-Newton 
was able to resolve some ~ 60 per cent of the 5-10 koV back- 
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I2OOII) . The Chandra Deep Fields North and South probe down 
to ~ 1.4-2.8 X 10"^*^ erg cm~^ s"^ (over 2-8 keV) and resolve 
~ 70-90 per cent of the background in this broad, hard band 
( lGiacconi etai]|200A iMoretti et al]|200l [Alexander et"^l2003t 
and references therein). 

The large resolved fractions found in the deep surveys appear 
to solve the problem of the origin of the XRB in terms of known 
sources; an extrapolation of the 0.5-2 keV band log A'^-log S dis- 
tribution points to an integrated flux consistent with the full XRB 
level. The substantial resolved fraction in the broad 2-10 keV band 
is not entirely consistent with this picture; a similar extrapolation of 
the 2-10 keV log A'^-log S distribution is, at most, only able to ac- 
count for some 93 per cent of the total XRB and is only marginally 
consistent with complete resolution (Moretti et al. 2003). Th e mos t 
recent XRB intensity measurement bv De Luca & Molendil |20o3) 
would make the fraction only ~ 80 per cent. 

To investigate i n detail the reso lved XRB fraction as a func- 
tion of energy, Worslev et alJ i2004) carried out a source-stacking 
analysis of ~ 700 ks of accumulated XMM-A^ewton exposure in the 
Lockman Hole. Source photometry was used to determine the re- 
solved fraction in a number of narrow energy bands from 0.2 to 
12keV. The authors found the resolved fraction to be > 80 per 
cent below 2 keV but that this decreased significantly at higher en- 
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ergies, falling to only ~ 50 per cent above ~ 8 keV. The failure to 
account for XRB in the harder bands, along with recent indications 
of steepening of X-ray source number counts at low fluxes, sug- 
gest that there may be an as-yet undetected population of highly 



recent synthesis models (e.g. 


Gilli et al. 


2001; 


Franceschini et al. 


bOOllGandhi & Fabiarj2003l:IUeda et ali2003: 


Comastri 2004). 



Although Chandra has low effective area at energies exceed- 
ing 7 keV, the Chandra Deep Fields (CDFs) probe to fainter flux 
limits than the XMM-Newton Lockman Hole observation. We use 
the CDFs to assess whether the contribution of the fainter sources 
found by Chandra are able to resolve the discrepancy and account 
for the missing background fraction. Using a similar photometric 
analysis we sum the flux from resolved CDF sources in a num- 
ber of narrow energy bands to determine the resolved fraction of 
the X-ray background. We also pres ent a re-analysis of th e XMM- 
Newton Lockman Hole observation iWorslev et alj|2004l) . carried 
out in identical energy bands for direct comparison with the CDFs. 
In both cases we attempt to quantify and correct for field-to-field 
variations and the missing flux which would arise from the very 
bright sources which are not sampled in these deep, pencil-beam 
surveys. 



2 SOURCE DETECTION AND PHOTOMETRY 
2.1 Chandra Deep Field Observations 

Chandra has performed the two deepest ever surveys of the 
X-r ay universe - the 2 M s Chandra Deep Field North (CDF- 
N; lAlexander e taP '2003") and the 1 Ms Chandra Deep Field 
South (CDF-S;^Giacconietal. 2002). The CDF-N probes down 
to ~ 2.5 X 10"^^ and - 1.4 x 10"" erg cm"^ s~^ in the 0.5-2 
and 2-8 keV bands respectively. Although the CDF-S only reaches 
down to equivalent limits (at S/N = 3) of ~ 5.2 x 10~^^ and 
~ 2.8 X 10~^®ergcm~^ s~^, these are both well beyond the 
10"^** er g cm~^ s~^ break in the logA'^-logS' distribution (see 
e.g.lJiasinger et al. 1998; Camoana et al. 2001). 

For both the CDF-N and CDF-S we use th e source cata- 
logues and photometry from lAlexander et alJi2003[) . who describe 
their approach thoroughly; we include a brief summary here for 
completeness. The 2 Ms CDF-N main point-source catalogue con- 
tains a total of 503 sources and covers a total solid angle of 
447.8 arcmin^. Conservative estimates place the number of falsely 
detected sources to be as high as 20-30 but the true number is prob- 
ably considerably lower (see section 2.3 m Alexander et al. 2003;. 
The CDF-S has 326 sources from a total field of 391.3 arcmin^. 

The local background estimate for each source was deter- 
mined in an annulus surrounding the source-extraction region. The 
mean number of background counts was computed assuming a 
Poisson distribution, as ni/no; where n\ is the number of pix- 
els containing a single count and no is the number containing no 
counts. For sources with less than 10'^ counts in the full 0.5-8 keV 
band, source-extraction apertures containing 90-95 per cent of the 
encircled energy fraction (EEF) were chosen (with the extracted 
counts corrected up to 100 per cent). For sources brighter than 
10'^ counts, extraction apertures of twice this size were chosen; the 
EEF in these cases is close to 100 per cent. The EEF corrections 
used to correct from aperture counts to total counts were calculated 
from the Chandra point spread function (PSF) dete rmined using the 
Chan dra X-ray Centre MKPSF software (refer to [Alexander et al] 
l2003h. 



For all sources with a sufficient number of counts, the ra- 
tio of counts between the 2-8 keV and 0.5-2 keV bands was 
used to determine an effective spectral index for a power-law 
model, allowing for Galactic absorption of (1.3 ± 0.4) x 10^" and 
(8.8 ± 4.0) X 10^'-^ cm"^ for the CDF-N and CDF-S respectively 
(Lockman 2003; Stark et al. 1992). For sources with an insufficient 
number of counts, the photon index is fixed at F = 1.4 (the ap- 
proximate average for all the sources and the slope of the XRB). 
The flux for each source was calculated in a number of photometric 
bands, regardless of whether or not the source is actually detected 
in each of the bands. In each case the counts-to-flux conversion 
is appropriate to the spectral index of the source. The photometric 
bands used were 0.5-1, 1-2, 2-4, 4-6 and 6-8 keV. 

Prior to our stacking analysis (described in sect ion 13. 1>. we 
applied the systematic flux correction determined bv lBauer et alJ 
(.2004.) for both the CDFs (see section 3 and Fig. 2 of their paper). 
These take into account Eddington bias as well as some further 
aperture/photometry effects that were not considered originally. 
The correction is negligible except in the harder (> 2keV) bands 
which each see an increase in total resolved flux of ~ 2-5 per cent. 

2.2 XMM-Newton Lockman Hole Observations 

XMM-Newton, whilst unable to match the X-ray resolution and ul- 
timate sensitivity of Chandra, has the critical advantage of sensi- 
tivity in the 8-12 keV regime. It was this capability, put into effect 
with ^ 700 ks of deep observatio ns in t he Lockman Hole (XMM- 
LH), that prompted lWorslev et ani2004h to investigate the resolved 
fraction of the XRB with energy. 

The 17 different individual observations that comprise the 
XMM-LH were combined and source detection carried out using 
sliding-box detection followed by maximum-likelihood PSF fit- 
ting. Visual inspection of the source candidates was used to ex- 
clude spurious detections. The analysis was restricted to the central 
314.2 arcmin'^ region of the field with the highest exposure time 
and lowest background level. 

In a similar approach to that taken in the CDFs, simple aper- 
ture photometry was performed on the exposure-corrected images 
in a number of narrow energy bands. A more detailed description 
of the following is given by Worslev et al. ( 2004). For each source, 
the total count-rates were extracted in both a circular aperture 
and a surrounding annulus. These were then used to calculate the 
background-corrected source count rate, taking into account cor- 
rections for the EEF. These corrections were calculated by integrat- 
ing an analytical model of the XMM-Newton PSF given in XMM- 
N ewton Sci ence Operations Centre calibration documents (refer 
to lWorsle'v'et al.2004) . Count-rate extraction apertures/annuli were 
chosen to the maximise the signal-to-noise ratio. 

Count-rate-to-flux conversion factors were computed assum- 
ing a fixed spectral slope of F = 1.4 and Galactic absorption with 
a column density of 5 x 10^^ cm~^ jLockman et alll986ib . This is 
not as rigorous as the approach taken in the CDFs where counts-to- 
flux conversion factors (ECFs) appropriate to the spectral slope of 
the source were used, but the larger count rate errors for the XMM- 
Newton sources prohibit this, and the assumption of a fixed F = 1.4 
is more reliable. Furthermore, the dependence of the ECFs on the 
spectral shape of a source is only weak because of the use of nar- 
row energy bands. The ECFs of all but the 0.2-0.5 keV band vary 
by < 5 per cent over the range F = 0.9-1.9 (the 0.2-0.5 keV band 
has a 15 per cent variation). The resultant uncertainty in the re- 
solved total flux from the source population is small in comparison 
with other errors. 
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Table 1. Parameters used in the XMM-LH photometric analysis. For each energy band the cut-out radii are quoted: ri is the radius of the source count-rate 
extraction aperture; r2 and are the inner and out radii of the background count-rate extraction annulus; refer to sections 2.3-2.4 and Fig. 1 in Worsley et al. 
(2004) for further details. Also quoted are the count-rate to flux energy conversion factors calculated for each band following the method described in section 
2.5 of that paper. 



Energy band Cut-out radii (arcsec) Weighted energy conversion factors 



(keV) 
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MOS 




(10 


count cm 


2 erg 1) 




ri 


r2 


'"3 


ri 


'"2 




PN 


MOS-1 


MOS-2 


0.2-0.5 


14 


30 


40 


12 


25 


35 


7.052 


1.263 


1.242 


0.5-1 


15 


30 


40 


13 


25 


35 


7.375 


1.757 


1.735 


1-2 


14 


30 


40 


12 


25 


35 


5.631 


1.994 


1.977 


2-4 


11 


30 


40 


10 


25 


35 


1.960 


0.7650 


0.7664 


4-6 


10 


25 


35 


8 


25 


35 


1.143 


0.4017 


0.4207 


6-8 


8 


20 


30 


6 


20 


30 


0.6694 


0.1347 


0.1411 


8-12 


6 


15 


25 


5 


15 


25 


0.1908 


0.02059 


0.02220 



In order to allow easy comparison with the Chandra data, the 
XMM-LH analysis was repeated using the same energy bands as in 
the CDFs with additional bands covering 0.2-0.5 and 8-12 keV. 
The re-processing of the data took advantage of the latest releases 
of the XMM-Newton SCIENCE ANALYSIS SYSTEM (SAS) V6.0.0 
and the X-ray spectral analysis package XSPEC V l 1.3.1. Addi- 
tion ally, rather than use the original source list (from lWorslev et all 
of 126 objects, we took advantage of the complete - and 
more rigorously compiled - Lockman Hole source catalogue of 
Mateos et al. (in prep.), which contains 156 sources in the same 
314.2 arcmin^ region. TableQgives the values of various impor- 
tant parameters used in the re-processing of the data. 
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3 RESULTS 

3.1 Resolving the X-ray background 

The total resolved flux in each energy band was calculated by sum- 
ming the measured flux for each source. Every source is included 
regardless of whether or not it is detected in the band. This was 
done for the CDF-N, CDF-S and each of the three XMM-Newton 
instruments in the XMM-LH. Appropriate corrections for Galac- 
tic absorption were made. These correction factors are 15 per cent 
in the XMM-LH 0.2-0.5 keV band; 2.6, 7.2 and 4.8 per cent in 
the 0.5-1 keV band for the XMM-LH, CDF-N and CDF-S respec- 
tively; and < 1 per cent for all of the harder bands. 

In order to calculate the resolved intensity, the flux from each 
source must be divided by a solid angle on the sky. In the case of 
the XMM-LH the sensitivity is approximately constant over the re- 
gion considered here and we simply take the solid angle to be a 
constant 314.1 arcmin^. The CDFs, however, show a substantial 
increase in sensitivity towards the centre of the fields. The total 
solid angles of 447.8 and 391.3 arcmin^ are only applicable to the 
brightest sources, whereas fainter sources are only detectable over a 
fr action of the field and th e actual solid angles are smaller. Fig. 19 
in [Alexander et alJ i2003h shows the solid angle as a function of 
0.5-8 keV source flux in the CDF-N. We use this function (and 
its equivalent for the CDF-S) to provide the appropriate solid angle 
when converting each source flux to an intensity. We impose a min- 
imum effective area of 10 per cent of the total to avoid catastrophi- 
cally magnifying the errors in the fluxes from the faintest sources. 

Fig. 0a) shows the resolved XRB fraction obtained in each 
energy band. The fractions are calculated using the most re- 
cent estimates of the total extragalactic XR B spectrum. For the 
1-8 keV band this is the power-law model of iDe Luca & Molendil 
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Figure 1. The fraction of the total extragalactic XRB intensity resolved 
by simply summing the fluxes of sources in the CDFs and XMM-LH. 
CDF-N and CDF-S data are shown as the black and grey crosses respec- 
tively; XMM-LH PN camera data are shown as dashed diamonds and both 
the MOS-1 and MOS-2 camera data are both shown using dotted dia- 
monds. The thi'ee panels show; (a) the result of summing all sources, (b) 
the result of excluding sources brighter than 5 X 10~^^ erg cm~^ in 
0.5-8 koV flux, and (c) the result of additionally excluding sources fainter 
than lO^^'' erg cm~^ in 0.5-8 keV flux. 
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i2004) . with a spectral slope of F = 1.41 and IkeV normalisa- 
tion of 11.6 keV cm"^ sr~^ keV"^ Above 8keV the back- 
ground spectrum starts t o turn over s l ightly and so here we use 
the analytical model of iGruber et alJ il999h - which takes the 
downturn into account - but renormalised to the same 2-8 keV 
XRB flux as found by De Luca & Molendi (2004). The spec- 
trum of the purely extragalactic background is not well con- 
strained below 1 keV since it is difficult to separate from the 
bright, diffuse, Galactic component. There is evidence to suggest 
that the spectrum o f the background steepens below ~ 1 keV; 
iRoberts & Warwick ^2001) have combined the results of several 
X-ray shadowing measurements and report a background intensity 
of 26.6 ± 4.8 keV cm"^ s"^ sr~^ keV~^ at 0.25 keV . A power- 
law fit between this value and the lPe Luca & Molendil m- 
tensity at 1 keV gives a spectral slope of F ~ 1.5-1.7 and we use 
this to estimate the background levels in the 0.2-0.5 and 0.5-1 keV 
bands. 

The stacking of source fluxes reveals a decline in the re- 
solved XRB fraction as a function of energy. Below ~ 2 keV 
the data diverge somewhat, presumably due to field-to-field vari- 
ations. Such variations arise through the cosmic variance in the 
actual XRB normalisation in different regions of the sky, along 
with the statistical noise which is a consequence of sampling a 
discrete population of sources in a small solid angle. Fig. 0b) 
shows the result of excluding sources with a 0.5-8 keV band flux 
> 5 X 10~^* erg cm~^ from the sum. While only a small 
number of sources have been removed (from 3-7 out of 156-503) 
the variation at the softest energies has been reduced considerably. 
This demonstrates that a large amount of the variation < 2keV 
is due to soft, bright sources. The XMM-LH data now show a 
slight decline in the resolved fraction whilst the CDFs show a slight 
rise, although both CDFs show significant drops in the 6-8 keV 
band. This difference between the XMM-Newton and Chandra frac- 
tions in the 4-6 keV band is due to sensitivity - the CDFs probe 
down to fainter fluxes where the typical source spectrum is harder 
than F = 1.4 (see e.g. [Alexander et alJl2003l: IStreblvanska et alJ 
l2003h . If the source lists are further restricted to remove the faint 
sources seen only by Chandra, with a minimum flux requirement 
of 10~^^ erg cva"^ s~^, then the spectral shapes of the sources are 
similar for both XMM-LH and CDF, see Fig.Qc). Of course, there 
will always remain differences between data-sets - Bar cons et alJ 
(2000) estimate the cosmic variation in XRB normalisation to be 
of order ~ 10 per cent for fields of < 1 deg^. 

3.2 Correcting for the bright-end population 

The brightest sources in the catalogues have 0.5-8 keV fluxes up to 
~ (1-3) X lO"^'^ erg cm^^ s^^. There is a non-negligible contri- 
bution to the whole-sky background by sources brighter than this; 
in fact, around 15-20 per cent of the total XRB is due to bright, rare 
sources that are not sampled in pencil-beam surveys. A 'bright-end 
correction' is necessary to recover this missing intensity. 

In order to remove some of the bright-end variation we trun- 
cate the source lists above 5 x 10~^^ergcm~^ (0.5-8 keV 
band). This is the same truncation as used in Fig.Qb) and corre- 
sponds to removing sources that only occur with densities of ~ a 
few in the ~ 0.1 deg^ CDF and XMM-LH fields. With the bright- 
end sources removed the expected missing flux can be estimated 
by integrating known log Ai'-log S distributions. These are avail- 
able in the broad 0.5-2 and 2-10 keV bands but not in the nar- 
row bands used in this work. In order to calculate the contribution 
in a narrow band it is necessary to assume the spectral shape of 



the sources. Here we assume a power-law spectrum but the typi- 
cal spectral slope of such a power-law can be a strong function of 
source flux. 

A number of studies point to sig nificant evolution of the 
spectral shape with flux lUedaetal. 1999; Fiore et all l2003t 
[Alexander et al.* 2003; Streb lyanska et al. 20()3), with sources 
brighter than '--^ 10~^* erg cm~^ having typical spectral 
slopes much softer than the F = 1.4 of the total background. 
Shallow surveys (Ue daetalj[l999l see e.g.) find that the bright- 
est sources, with X-ray fluxes > 10~^'^ erg cm~^ s~^, h ave spec- 
tral slopes of F ~ 1.8-2.2. Istreblvanska et alJ J2003h fitted a 
power-law spectrum to the stacked spectra of sources of a 500 ks 
XMM-Newton observation in the CDF-S. Fitting the hard band 
(1-8 keV) with a power-law reveals a strong correlation between 
spectral slope and flux, with F ~ 1.8 for ^ 10~^^ erg cm~^ 
(2-10 keV flux) decreasing to F ~ 1.5 at ^ 10"^"^ erg cni~^ 
and F ~ 1.2 at ~ 10"^^ erg cm"^ s"\ [Alexander etall <2003h 
see sources with F ~ 1 at fluxes < 10~^® erg cm~^ in the 
CDF-N. 

We haven taken a linear fit to the evolution of spec- 
tral index with flux bas ed on the relationship as observed by 
Istreblvanska et all i2003h . although we conservatively impose a 
maximum value of F = 2. Such a model is only an approximation 
to the actual dependence of spectral shape with flux. It does, ho w- 
ever, follow the trends observed by Streblyans ka et alJ i2003h in 
the CDF-S (using both XMM-Newton and Chandra observations) 
and Ueda et al. ( 1999) using the ASCA Medium-Sensitivity Sur- 
vey, and whilst crude, it is important to take into account the fact 
that the bright-end correction is due to sources which are signifi- 
cantly softer than those which we have already stacked. Given our 
model, the bright-source contribution in each of the 2-4, 4-6, 6-8 
and 8-12 keV narrow bands can then be calculated, using an ap- 
propriate spectral index at each source flux, from the 2- 1 keV 
flux contribution obtained by integrating the iMoretti et alJ i2003h 
hard band log A'^-log S distribution (the most complete published 
to date). 

Correcting the soft bands (0.2-0.5, 0.5-1 and 1-2 keV) is 
more complex, since the shape of the backg round is not well 
constrained. The source stacking approach of IStreblvanska et alJ 
(200a|) in the 0.4-8 keV range revealed a systematically softer 
spectrum than that seen for 1-8 keV (an increase in F of ~ 0.1). 
No studies cover the 0.5-2 keV range ideally required for the soft 
band bright-end correction so we take the same linear fit to the 
evolution of spectral index as used in the hard band, but with an 
increase in F by 0.1. We use the Moretti et al. (2003) 0.5-2 keV 
soft band log N-\ o^ S distribution. Th ere is now some evidence 
to suggest that the iMoretti et'aP i2003h log A'^-log S slope is too 
steep at the bright end compared to resu l ts from the RIXOS , 
RBS and NEPS surveys ^ Mason et alJbOOOHSchwope et alJbOOOt 
^ioia et al. 2003); thus, our bright-end correction in the soft band 
could potentially be an underestimate. 

Fig. |2| shows the bright-end corrected resolved background 
fractions. The XMM-LH, CDF-S and CDF-N are all consistent 
with a similar resolved fraction of ~ 70-90 per cent up to ~ 4 keV. 
The differences are of order ~ 10 per cent which is not inconsistent 
with what would be expected from cosmic variance between the 
fields (see lBarcons et al]l2000l) . XMM-Newton sees a drop in the 
resolved fraction in the 4-6 keV band whereas the CDFs are still 
consistent with a high resolved fraction of ~ 80-90 per cent. This 
difference is due to the faint, hard sources which are not detected in 
the XMM-LH; if these sources are removed from the Chandra data 
then all three data-sets show agreement in the 4-6 keV fraction. 
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Most importantly, the picture is different at high energies; XMM- 
Newton and both Chandra deep fields only resolve ~ 50-70 per 
cent of the background in the 6-8 keV band, much less than the 
fraction at lower energies. The 8-12 keV band, as probed by XMM- 
Newton, is only ~ 40-60 per cent resolved. 

We find the total resolved fluxes in the 0.5-2 
and 2-8 keV bands to be (6.9 ± 0.2) x 10"^^ and 
(1.5±0.1) X 10""ergcm-^ for the CDF-N. 

These compare well with the (6.7 ± 0.3) x 10"^^ and 
(1.66 ± 0.05) X 10~^^ erg cm~^ reported in iBauer et all 
yOOJ) from a detailed analysis of the combined CDF 
number counts. The fluxes can also be compared to the 
(7.1 ± 0.3) X 10"" and (1.43 ± 0.08) x 10"" erg cm~^ 
determ ined integrating the log A'^-log S curves of Mo retti et alj 
i2003 ^ (although note that these do not contain the 2 Ms extension 
to the original 1 Ms CDF-N). 

We find the total 0.5-2 and 2-8 keV resolved fractions are 
~ 85 and ~ 80 per cent for the CDF-N (using the total XRB spec- 
trum we describe in section lTTl . Such high fractions are often used 
to claim that the XRB is almost completely resolved into discrete 
sources even in the broad 2-10 keV band. Our results suggest that 
stating a resolved fraction in this manner is misleading, since it is 
the high resolved fractions in the 2-6 keV range which are domi- 
nating the total, whilst the resolved fraction beyond 6 keV is sig- 
nificantly lower. 



4 DISCUSSION 

Three independent observations now point to a significant reduc- 
tion in the resolved fraction of the X-ray background at the hardest 
energies. XMM-Newton starts to see a downturn above ~ 4 keV 
whilst Chandra is sensitive to the fainter, harder sources, which 
fill in some of the missing fraction and push the downturn to en- 
ergies > 6 keV. The missing fraction has the spectral signature of 
highly obscured AGN and could be evidence of the large popula- 
tion of heavily absorbed objects which are t hought to account for 
the 3 keV peak of XRB intensity (see e.g. iMaioUno et alJll998l 
l2003h . 

To quantify the possible nature of an undetected popula- 
tion, we modelled the spectral shape of the unresolved back- 
ground as a function of redshift and intrinsic absorption col- 
umn density. A grid of spectra was built in the range z = 0.1-3 
and Nh = 10^^-10^^ cm~^. XSPEC Vll.3.1 was used to gen- 
erate the spectra using a PEXRAV model I Maedziarz & Zdziarski 
1 19951) . with a F = 2 power-law and a n R = 1 reflection com- 

S ionent (see e.g. iMa liziaetalJ l2003h . multiplied by ZWABS 
Morrison & McCammorT^SsT photoelectric absorption at the 
source redshift. 

We computed the goodness-of-fit between each model spec- 
trum and the residual XRB spectrum (i.e. the difference between 
the resolved and total backgrounds). Fig. |3| shows the confidence 
contours in the redshift and absorption column density required 
to account for the shape of the missing AGN population for the 
XMM-LH. This is a rather simplistic approximation to the spec- 
tra of highly obscured sources; specifically, we neglect the effects 
of iron K emission and scattered-flux components, as well as as- 
suming that the missing sources occur at a fixed redshift and with 
a fixed absorption column density. Nevertheless, some generalised 
conclusions about the nature of the unresolved population can be 
drawn. 

The XMM-LH, CDF-N and CDF-S all show similar contour 
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Figure 2. The upper panel shows the total extragalactic XRB (thick black 
crosses) along with the intensity resolved into sources in the CDF-N 
(black crosses), CDF-S (grey crosses) and combined PN/MOS-l/MOS-2 
results in the XMM-LH (black diamonds). The resolved intensities have 
been corrected for the bright-end population by removing all sources with 
0.5-8 koV flux > 5 x 10"^'' erg cm~^ s^^ and then adding the missing 
flux from sources brighter than this by integrating log Af-log S distribu- 
tions (refer to section IXll . The lower panel shows the resolved intensities 
as fractions of the total extragalactic XRB in each energy band. 



diagrams for the characteristics of the missing source population. 
The XMM-LH data, which extends to the 8-12 keV band, shows 
the most constrained set of contours although the degeneracy be- 
tween 1ow-z/1ow-7Vh and high-z/high-A'^H is clear. The best-fitting 
point is at 2 ~ 0.8 and Nu ~ 4.5 x 10^^ cm"^. For the CDF-N 
and CDF-S the lack of data above 8 keV results in much more re- 
laxed contours although there is again some preference for sources 
at z ~ 0.5-1.5 with A^h ~ 10^^-10^" cm"^ Figure IH shows the 
total and resolved levels of the XRB in vli, as well as the intensity 
that would be contributed by the best-fitting model of an obscured 
population. 

We conclude that whilst the data are unable to constrain tightly 
the spectral shape of the missing XRB sources, the most plausi- 
ble population would seem be of objects at redshifts of ^ 0.5-1.5, 
with unabsorbed luminosities of < 5 x lO'*'^ erg s~^, but heavily 
obscured by column densities of ~ 10^^-10^* cm"'^. The CDFs 
probe deeper than the XMM-LH but still fail to detect the missing 
sources in the 2-4 keV band, or our stacking analysis would have 
revealed the higher energy emission. Any soft starburst or scattered 
emission in these objects must be small and this may suggest a high 
covering fraction for the absorption. 
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Figure 3. Contour plot of fits to the spectral shape of the missing XRB 
spectrum for the XMM-LH data in the 2-12 keV range. 68, 90 and 95 
per cent confidence contours are shown. The best-fitting point is indicated 
at 2 ~ 0.8 and A'^h ~ 4.5 X 10^^ cm^^. The grey contours indicate the 
maximum, unobscured, rest-frame 2-10 kcV luminosity (ergs^^) of the 
source population required such that the sources remain below the sensitiv- 
ity limit (estimated at ~ 1.5 X 10~^^ erg cm~^ in a 4-8 keV band). 

In order to provide the missing fraction of the hard X-ray 
background these sources must have a density on the sky of 
> 2800 deg^^. This would correspond to a lower limit of around 
350 undetected AGN in the CDF-N. There must therefore be > 3 
times more undetected, obscured sources than detected, unobscured 
ones, in the 4-8 keV band. This requires that the sources in the new 
population have a high covering fraction of obscuring matter. 

From the absorption-corrected flux we can estimate the 
black hole mass density of this obscured population to be 
(6-9 X 10'*) eZ\ «_i Mq Mpc~^; where e_i is the accretion ef- 
ficiency in units of 0.1 and k_i is a bolometric correction in units 
of 0. 1. We have used a similar calculation to Fab ian & Iwasawal 
il99ft') although note t hat the b olometric correction is ~ 10-15 for 
our lower-luminosity jFabiaijE0 04y The missing population thus 
accounts for ^ 10-20 per cent of the local black hole mass density 
of 4-5 X lO'^ Mq Mpc"^ (see e.g. Shankar et al. 2004), the rest of 
which is due to bright AGN and, in particular, quasars, where the 
bolometric correction is ~ 30 or mor e, and where a typical reds hift 
of 2 ~ 2 is more appropriate (refer to lFabian & Iwasawall999l and 
references therein). 

The steepness of the 5-10 keV source number counts, with no 
evidence for a break to a flatter slope at the faint end, also indicates 
that there remains a substantial number of obje cts which have yet 
to be found in th is band iHasinaer et al. 2001 : iRosati et aljl2002t 
iBaldi et nl 'lQQ't). Recent optical and infra-red data from the HST 
and Spitzer telescopes, acquired as part of the GOODS project, 
are also prov iding evidence for l arge numbers of obscured AGN 
at 2 > 1 (e.g. lTreister et all2004b . Deep, multi-wavelength efforts 
such as GOODS may well be the best way to identify the missing 
population; confirmation of Compton-thick objects will ultimately 
require X-ray missions with high sensitivities beyond 10 keV, such 
as the next generation of focusing hard X-ray telescopes. 



5 CONCLUSIONS 

• Whilst the XRB is ~ 85 and ~ 80 per cent resolved in the 
broad 0.5-2 and 2-10 keV bands respectively (depending on XRB 
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Figure 4. The upper curve shows the total extragalactic XRB (as de- 
scribed in section ITT] in ulv The intensity resolved into sources in the 
CDF-N, CDF-S and XMM-LH are shown as black crosses, grey crosses 
and black diamonds respectively. The lower curve shows the spectrum 
of the best-fitting model for the missing AGN population with z = 0.8 
and Nn = 4.5 X 10^^ cm~^. Sharp features (such as that due to iron at 
4 keV), is a consequence of our assumption that the missing population is 
composed entirely of AGN with this spectrum. In reality the missing pop- 
ulation will be mixture of sources at different redshifts and with different 
levels of absorption; any sharp features will be smeared-out. 



normalisation), it is only ~ 60 per cent resolved above ~ 6 keV 
and ^ 50 per cent resolved above 8 keV. This decrease in re- 
solved fraction is seen in both the CDFs and the XMM-LH. 

• The decrease in resolved fraction as seen in XMM-LH occurs 
at energies > 4 keV whilst in the CDFs the drop occurs at > 6 keV. 
This difference is due to faint, hard sources which are detected in 
CDFs but not the XMM-LH. All three surveys are unable to ac- 
count for almost half of the X-ray background above ~ 6 keV. 

• The missing fraction has a spectral shape that is consistent 
with that which would be expected from a population of faint, heav- 
ily obscured AGN located at a redshift of ~ 0.5-1.5 and with in- 
trinsic absorption columns of ~ 10^^-10^* cm^"^. 

• The deep CDF data show that any such obscured population 
is not detected in the 2-4 keV band, or our stacking analysis would 
have revealed the higher energy emission; consequently, any soft 
scattered or starburst emission from these objects must be small. 
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